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ABSTRACT 

The spectral index of synchrotron emission is an important parameter in understanding 
the properties of cosmic ray electrons (CREs) and the interstellar medium (ISM). We deter¬ 
mine the synchrotron spectral index (ant) of four nearby star-forming galaxies, namely NGC 
4736, NGC 5055, NGC 5236 and NGC 6946 at sub-kpc linear scales. The a^ was determined 
between 0.33 and 1.4 GHz for all the galaxies. We hnd the spectral index to be flatter (> -0.7) 
in regions with total neutral (atomic -H molecular) gas surface density, Sgas > 50 Mopc“^, 
typically in the arms and inner parts of the galaxies. In regions with Sgas < 50 Mopc^^, es¬ 
pecially in the interarm and outer regions of the galaxies, the spectral index steepens sharply 
to < -1.0. The flattening of ant is unlikely to be caused due to thermal free-free absorption 
at 0.33 GHz. Our result is consistent with the scenario where the CREs emitting at frequen¬ 
cies below ~ 0.3 GHz are dominated by bremsstrahlung and/or ionization losses. For denser 
medium (Sgas > 200 Mopc^^), having strong magnetic fields (~ 30 juG), ant is seen to be flat¬ 
ter than -0.5, perhaps caused due to ionization losses. We And that, due to the clumpy nature 
of the ISM, such dense regions cover only a small fraction of the galaxy (< 5 percent). Thus, 
the galaxy-integrated spectrum may not show indication of such loss mechanisms and remain 
a power-law over a wide range of radio frequencies (between ~ 0.1 to 10 GHz). 

Key words: galaxies : ISM - (ISM:) cosmic rays - ISM : magnetic fields - galaxies: ISM - 
radiation mechanisms : non-thermal - radio continuum : galaxies 


1 INTRODUCTION 

The synchrotron spectral index of the cosmic ray electrons (CREs) 
in normal star-forming galaxies is an important parameter in un¬ 
derstanding their energy loss/gain mechanisms in the interstellar 
medium (ISM). Studying spatially resolved properties of the syn¬ 
chrotron (also referred to as non-thermal) spectral inde)|3 (ttnt) 
can throw meaningful insights into the properties of the ISM and 
thereby understand the nature of multi waveband emission from 
galaxies. 

It is believed that the CREs are produced by diffusive shock 
acceleration (DSA) of thermal electrons up to energies of few PeV 
in shock fronts of supernova remnants (Bell 2004; Jones 2011; Ed- 
mon et al. 2011; Kang, Edmon & Jones 2012). The CREs are in¬ 
jected into the ISM with a power-law energy spectrum n(E)dE oc 
where, n(E) is the number density of CREs with energies in 
the range E and E + dE. The energy spectral index, y, of CREs is 


* E-mail: abasu@mpifr-bonn.mpg.de 

* Spectral index, a, is defined as: Sy oc v®, where Sy is the radio continuum 
flux density at a frequency v. 


related to a„, as y = 2a„t - 1 and can be measured directly from 
observations. Depending on the CR acceleration mechanism and 
efficiency, y is estimated to lie in the range -2 and -2.4 (Biermann 
& Strom 1993; Bogdan & Volk 1983; Blandford & Eichler 1987; 
Bell 1978), which corresponds to in the range -0.5 and -0.7. 

The CREs that are injected into the ISM are subjected to en¬ 
ergy dependent losses and hence is expected to be modified 
smoothly over a frequency range. The main energy loss mech¬ 
anisms involving ionization losses, relativistic bremsstrahlung 
losses, synchrotron and inverse-Compton losses alter the radio fre¬ 
quency spectrum in a specific way. These loss mechanisms have 
different dependence on CRE energy (see Longair 2011) and thus 
they are important in different parts of the spectrum. Ionization and 
relativistic bremsstrahlung losses typically affect CREs having en¬ 
ergies < 1 GeV while synchrotron and inverse-Compton losses af¬ 
fect CREs having energies > 2 GeV. These energies correspond to 
~ 0.2 GHz and 1 GHz, respectively, for typical magnetic fields 
of ~ 10 yUG. 

At higher radio frequencies, synchrotron and inverse- 
Compton losses lead to a cut-off in the spectrum. However, the form 
of this cut-off depends on the pitch angle between the magnetic 
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field and CREs and the process of particle injection. In the case 
of steady continuous injection (Cl) of CREs, the spectrum steep¬ 
ens by 0.5, i.e., (Aa = -0.5) (Pacholczyk 1970). Eor injection of 
particles at a single epoch and considering constant pitch angles of 
individual electrons, according to the Kardashev-Pacholczyk model 
(KP; Kardashev 1962; Pacholczyk 1970) the spectrum falls off as 
a power-law with index 4ai„j/3 - 1, i.e., Aa = Q-inj/3 - 1. Here, 
Qfinj is the initial power-law index of the particles. The Jaffe-Perola 
model (JP; Jaffe & Perola 1973) assumes rapid isotropization of 
the pitch angle distribution leading to an exponential cut-off in the 
CRE spectrum. Eurther, energy dependent diffusion losses can also 
steepen the spectrum by 0.25, i.e., Ao- = -0.25 (Condon 1992). 
However, such losses would affect the CREs having high energies 
(> 10 GeV). 

On the other hand, at low radio frequencies, the ionization 
losses have the effect of flattening a„f Under steady continuous 
injection of electrons, the equilibrium spectra is flatter by 0.5, i.e., 
Ao-n, = -1-0.5. Eor single injection events, the resulting spectrum 
will depend on the ratio of timescales of injection and ionization 
loss. Bremsstrahlung losses mildly affects the CRE spectrum and 
therefore does not alter o-n, significantly (Attn, =» 0) (see Longair 
2011), both in the case of continuous and single injection. Thus 
comparing the value of the observed with respect to the typi¬ 
cal injection values of-0.5 to -0.7 one can broadly get an idea 

of the dominant CRE energy loss mechanism. Eurther, due to ther¬ 
mal free-free absorption, the synchrotron spectrum is expected to 
turn over at even lower frequencies below ~ 0.1 GHz. These have 
a combined effect of modifying the initial power-law spectrum into 
a curved spectrum with flatter slope at low frequencies (< 1 GHz) 
and steeper slope at high frequencies (> 1 GHz). 

The radio continuum emission from galaxies broadly origi¬ 
nates due to two main emission mechanisms - 1) synchrotron emis¬ 
sion and 2) thermal bremsstrahlung emission. Synchrotron emis¬ 
sion contributes more than 90 percent to the total radio emission 
at low frequencies (< 1 GHz) (Niklas, Klein & Wielebinski 1997; 
Basu et al. 2012). Contribution from thermal emission to the to¬ 
tal radio emission is significant at higher frequencies. Eor example, 
at 10 GHz, ~ 40 percent of the total emission is thermal in ori¬ 
gin (Gioia, Gregorini & Klein 1982). The thermal emission is also 
characterized by a power-law spectrum having thermal spectral in¬ 
dex ttth = “0.1 and is significantly flatter than that of the The 
composition of these two types of emission gives rise to an overall 
power-law spectrum (a) which is flatter than ornt- The presence of 
thermal emission makes it difficult to interpret energy spectra of the 
CREs and therefore it is necessary to subtract its contribution from 
the total emission. Note that, in this paper we distinguish a and 
as the spectral index of the total radio and synchrotron emission, 
respectively. 

A host of radio continuum surveys has been undertaken in the 
past to study the broadband integrated spectra of galaxies and un¬ 
derstand the various energy loss mechanisms. Surprisingly, most 
of the studies showed a power-law spectrum over broad range of 
frequencies. The power-law slope is found to be similar, having a 
value-0.8 with narrow dispersion, typically less than 20 per¬ 

cent. The power-law index, a between 1.4 and 5 GHz was found to 
be -0.85 (Sramek 1975); between 0.4 and 10.7 GHz, a was found 
to be -0.71 (Klein & Emerson 1981) and -0.74 (Gioia et al. 1982). 
Based on broadband fitting of the radio continuum spectrum be¬ 
tween 0.6 to 10 GHz for 41 spiral galaxies, Duric, Bourneuf & 

Gregory (1988) found to peak at-0.8. In a similar study of 

74 Shapely-Ames galaxies between 0.4 and 10.7 GHz, Niklas et al. 
(1997) found the mean value of am to be -0.83. In fact, for a much 


wider range of radio frequencies, between 26.3 MHz and 22.8 GHz, 
Mulcahy et al. (2014) did not find any evidence of curvature in the 
radio continuum spectrum for the galaxy M51. The spectrum is 
consistent with a single power-law with a = -0.79 ± 0.02. 

Recently, Marvil, Owen & Eilek (2014), in a statistical study 
of 250 galaxies, found evidence of curved radio continuum spectra 
between 74 MHz and 4.85 GHz. However, there is no strong ev¬ 
idence of spectral curvature when studied for integrated flux den¬ 
sities of individual galaxies except for a few cases, such as M82 
(Condon 1992; Adebahr et al. 2013), NGC 4631 (Pohl & Schlick- 
eiser 1990), NGC 3627 and NGC 7331 (Paladino, Murgia & Orrii 
2009), etc. At much lower frequencies of 57.5 MHz, Israel & Ma¬ 
honey (1990) found evidence of systematically lower flux densities 
when compared to extrapolated flux densities from higher frequen¬ 
cies assuming a single power-law. They concluded that smoothly 
distributed diffuse ionized gas throughout the galaxies is unlikely 
to cause thermal free-free absorption and requires more clumpy 
medium of well-mixed non-thermal emitting and thermally absorb¬ 
ing gas. 

Although the global spectral index, for both total and non- 
thermal emission, of galaxies does not vary and is consistent with 

ttm-0.8, locally they show large variations. The spectral index is 

observed to vary significantly radially, being flatter and close to the 
injection value (am ^ -0.5 to -0.6) towards the central regions and 

steepens towards the outer parts (a^,-0.9 to -1.2) (Basu et al. 

2012; Beck 2007; Tabatabaei et al. 2007; Paladino et al. 2009). The 
various CR energy loss mechanisms depend on local physical pa¬ 
rameters of the ISM such as the magnetic field strength, the number 
density of neutral gas, energy density of the photon field, shocked 
regions, etc. Moreover, the magnetic field strengths and gas density 
are believed to be coupled due to magnetohydrodynamic turbulence 
in the ISM (Chandrasekhar & Eermi 1953; Cho & Vishniac 2000; 
Groves et al. 2003). Eurther, injection rate of CR particles depends 
on the star formation rate which is related to the gas density through 
the Kennicutt-Schmidt law (Kennicutt 1998). Local ISM densities 
can therefore, directly or indirectly shape up the spectral index dif¬ 
ferently giving rise to large variations. It is therefore imperative to 
study the spatially resolved non-thermal spectral index and com¬ 
pare it to the local ISM gas densities. In this paper we qualitatively 
discuss the effects of gas surface densities on the non-thermal spec¬ 
tral index and the most likely scenario of CR energy loss prevalent 
in spatially resolved galaxies. 

The paper is organized as follows: In Section 2 we present 
the various source of the data used. The galaxies studied are NGC 
4736, NGC 5055, NGC 5236 and NGC 6946. We present our re¬ 
sults in Section 3 and discuss them in Section 4. The conclusions 
based on this work are summarized in Section 5. 


2 DATA AND ANALYSIS 

In this work we study four nearby galaxies, namely, NGC 4736, 
NGC 5055, NGC 5236 and NGC 6946. The galaxies are chosen 
from Basu et al. (2012) and were observed using the Giant Metre- 
wave Radio Telescope (GMRT) at 0.33 GHz. To measure the spec¬ 
tral index, we used archival data at a higher radio frequency near 1.4 
GHz. Table[T]lists the salient features of the sample and the sources 
of the archival data. As pointed out earlier, the radio continuum 
emission originates from both synchrotron emission and thermal 
free-free emission. Owing to the steep spectral index, a,,, < -0.5, 
the synchrotron emission is expected to dominate over the thermal 
emission, especially at lower (< 1 GHz) radio frequencies. How- 
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Table 1. The sample galaxies. 


Name 

Type 

Angular 

i 

Distance 

CO 

Hi 

Radio 



size (D 25 )(') 

n 

(Mpc) 



(1.4 GHz) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(V) 

(8) 

NGC 4736 

SAab 

11.2x9.1 

41 

4.66* 

HERACLES 

THINGS 

Westerbork* SINGS'* 

NGC 5055 

SAbc 

12.6x7.2 

59 

9.2+ 

HERACLES 

THINGS 

Westerbork SINGS'* 

NGC 5236 

SABc 

11.2x11 

24 

4.5L 

NRAO 12 m 

THINGS 

VLA“ CD array+ 

NGC 6946 

SABcd 

11.5x9.8 

33 

6.8+ 

HERACLES 

THINGS 

VLA C-l-D array* 


Column 3 lists the optical diameter measured at the 25 magnitude arcsec^^ contour. The inclination angles (i; 0° is face-on) are listed in Column 4. Distances 
in column 5 are taken from: * Karachentsev et al. (2003), ^ Karachentsev et al. (2002), ^ Karachentsev, Shaiina & Huchtmeier (2000) and the NED t. 
Columns 6 and 7 list the data used to trace the molecular and atomic gas, respectively. Column 8 lists the sources of archival data at 1.4 GHz: Braun et al. 

(2007), ^ VLA archival data using CD array (project code: AS325), ® VLA map obtained by combining archival data from C and D array (Beck 2007). 

“ The Very Large Array (VLA) is operated by the NRAO. The NRAO is a facility of the National Science Loundation operated under cooperative agreement 

by Associated Universities, Inc. 

* The Westerbork Synthesis Radio Telescope (WSRT) is operated by the Netherlands Loundation for Research in Astronomy (NLRA) with financial support 

from the Netherlands Organization for scientific research (NWO). 


Table 2. Resolution of the available maps in arcsec^. 


Name 

0.33-GHz 

1.4-GHz 

Hi 

CO 

^nt 

maps 

Linear 
scale (kpc) 

NGC 4736 

13x12 

19x12.5 

10.2x9.1 

13.4x13.4 

20x20 

0.45 

NGC 5055 

17x10 

18.5x12.5 

10.1x8.7 

13.4x13.4 

20x20 

0.90 

NGC 5236 

16x12 

26x14 

15.2x11.4 

55x55 

26x14 

0.55 

NGC 6946 

12x11 

15x15 

6x5.6 

13.4x13.4 

15x15 

0.50 


ever, due to high star-formation activity and/or certain bright Hit 
regions in galaxies, the contribution of the thermal emission can be 
significant locally (Basu et al. 2012). The contamination due to the 
thermal component of the radio emission, which has significantly 
flatter spectral index of -0.1, affects the synchrotron spectral index, 
making it flatter than actual by more than 30 percent (see Figure 6 
of Basu et al. 2012). This necessitates a thorough separation of the 
thermal emission from the total radio emission for robust determi¬ 
nation of ttnl- 

We used the recent method of estimating the thermal compo¬ 
nent of the radio emission given in Basu & Roy (2013). The ther¬ 
mal emission is determined by combining Hor and A24 /rm emis¬ 
sion from galaxies as its tracer. This method allows us to probe the 
non-thermal emission from galaxies at linear scales of 0.4-1 kpc, 
depending on the distance of the galaxies. Tablel^lists the angular 
resolution of the available maps. The overall thermal fractioijl 
at 0.33 GHz (/th.o.saoHz) was found to be less than 5 percent while 
at 1.4 GHz (/th,i.4GHz) it is found to be ~ 10 percent (Basu & Roy 
2013; Basu et al. 2012). However, locally, /th.o.ssoHz and /th,i.4GHz 
are as high as 15 percent and 30 percent, respectively. 

To estimate the surface mass density of total neutral gas (Egas) 
in the galaxies, we used Hi and CO line emission as tracers of 
the neutral atomic and molecular gas components, respectively. We 
used the natural-weighted Hi moment-0 maps from the Hi Nearby 
Galaxy Survey (THINGS; Walter et al. 2008) to trace the atomic 
gas for all our sample galaxies. The galaxies were observed with 
high angular resolution of ~ 10 arcsec (see Table |2j and velocity 
resolution < 5.2 km s“' using the VLA. We used CO(J2^1) emis- 


^ Thermal fraction (/ih.v) at a radio frequency y is defined as: /ih.v = 
5th,v/5tot,v Here, 5th,v and Stot.v are the estimated thermal emission and 
the total radio emission, respectively. 


Sion from the galaxies, NGC 4736, NGC 5055 and NGC 6946 as a 
tracer of the molecular gas in the ISM. We used the CO(I2—>1) 
moment-0 maps from the HERA CO Line Extragalactic Survey 
(HERACLES; Leroy et al. 2009). The galaxies were observed with 
an angular resolution of ~ 13 arcsec and velocity resolution of 
2.6 km s“* using the IRAM 30-m telescope. For NGC 5236, we 
used CO(Il—>0) line emission as the tracer of molecular gas ob¬ 
served with the NRAO 12-m telescope (Crosthwaite et al. 2002). 
The CO(Jl—>0) moment-0 map has an angular resolution of 55 arc¬ 
sec and velocity resolution of 5.2 km s“*. 

We computed total neutral gas surface mass density maps on a 
pixel-by-pixel basis using the surface mass density maps of atomic 
(Ehi) and molecular gas (Ehj) as, Egas = L36x (Ehi lEhj). The 1.36 
factor is to account for the presence of Helium. The Hi gas surface 
mass density was estimated using, 

Ehi(MoPC- 2) = 0.015 cos i /hi (K km s"'). (1) 

Here, i is the inclination of the galaxy, defined such that i = 0 is 
face-on and /hi is the line-integrated intensity. 

Similarly, the CO(I2—>1) line-integrated intensity, /co 2 ^ii of 
the HERACLES maps were converted to the molecular gas surface 
mass density using, 

EnjfMopc^^) = 5.5 cos i /coj^i (K km s“'). (2) 

This assumes a constant CO-to-H 2 conversion factor, Xco = 2 x 
10^°cm“^ (K km s“*)“* for the CO(Jl—>0) transition and a line ra¬ 
tio, /co 2 ^i/fcoi^o “ Tf*® line-integrated flux density, Scoi^o’ 
of the CO(Jl—>0) moment-0 NRAO 12-m map for the galaxy NGC 
5236 was converted to molecular gas mass, A/hj (Young et al. 
1989), 

M„,(Mo) = 1.1 X lO^D^p, cos i Sco.^ody km s-‘). (3) 

Here, Dmpc is the distance to the galaxy in Mpc. The molecular gas 
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NGC 4736 



NGC 5236 



Frequency (GFIz) 


NGC 5055 



NGC 6946 



Figure 1. Integrated radio continuum spectrum of the sample galaxies. The black dots show the total observed flux density and the red squares show the 
estimated non-thermal flux density. The solid and dashed lines show the power-law fit to the total and non-thermal emission respectively. The green dotted 
lines show the extrapolated thermal emission. 


mass was converted to Ehj by EhjCMqPC^^) = x 

Tpix (pc^)]. Here, Xpix and ypix are the pixel sizes of the map ex¬ 
pressed in pc at the distance of the galaxy. 


3 RESULTS 

Figure [T] shows the galaxy integrated radio continuum spectrum 
for the four galaxies between 0.33 and 10.7 GHz. The black dots 
represent the total observed flux densities obtained from vari¬ 
ous sources in the literature (see Table 3 of Basu et al. 2012). 
The higher frequency (> 5 GHz) data points were mostly mea¬ 
sured using single dish observations and hence, do not suffer from 
missing flux issues. The broad-band spectrum was fitted using 2 
types of models in log-log space; 1) a single power-law and 2) 
a second order polynomial that could account for any curvature 
in the spectrum, if present. The single power-law has the form, 
log S V = c + a log y. The second order polynomial has the form, 
log Sy = c -I- a logy + /3(logy)^. Here, c, a and are the fitting 
parameters. 


For all the galaxies, the best fit model is well represented by 
a single power-law shown as the solid black lines in Figure [T] The 
power-law index, a, is found to be -0.62 ± 0.05, -0.82 ± 0.01, 
-0.74 ± 0.02 and -0.70 ± 0.02 for the galaxies NGC 4736, NGC 
5055, NGC 5236 and NGC 6946, respectively. The red squares 
show the non-thermal flux densities estimated by subtracting the 
thermal emission (shown as the green dotted lines) at the respec¬ 
tive frequencies. The non-thermal emission spectrum for all the 
four galaxies can also be well fitted by a power-law having index 
steeper than the total emission. The non-thermal power-law index, 
ttni, is estimated to be -0.70 ± 0.05, -0.91 ± 0.01, -0.80 ± 0.03, 
-0.82 ± 0.02 for the galaxies NGC 4736, NGC 5055, NGC 5236 
and NGC 6946, respectively. Within the uncertainty of the mea¬ 
surements, we do not find any steepening of the overall spectrum 
towards higher frequencies (> 2 GHz) as expected due to syn¬ 
chrotron and/or inverse-Compton losses. Towards, the lower fre¬ 
quencies (< 0.5 GHz), we do not observe any flattening due to 
ionization losses and/or thermal free-free absorption. Our result is 
similar to what was observed for the galaxy M51 over a much wider 
range of frequencies (Mulcahy et al. 2014). 
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Figure 2. Non-thermal spectral index (am) versus the total neutral gas (atomic + molecular) surface density (Sg^s) in units of Mopc^^. The dotted lines shows 
am = —0.5. The blue stars shows the galaxy averaged values. We show representative error bars for some of the points covering the range of Sgas for each 
galaxy. 


In Figure 12 we plot the variation of the am with total surface 
density of the neutral (atomic+molecular) gas (Sgas) for all the four 
galaxies, am and Zgas is computed within one synthesized beam 
given in Table |2 To ensure independence, each point is measured, 
roughly separated by one synthesized beam. The am was computed 
using the non-thermal maps between 0.33 and 1.4 GHz consider¬ 
ing pixels having signal-to-noise ratio more than 30“ at both the 
frequencies. Each beam corresponds to 0.5 - 1 kpc in linear size 
depending on the distance of the galaxies (see Table |2j- For the 
galaxy NGC 5236, due to poor resolution of the available CO map, 
the linear resolution was limited to 1.2 kpc. 

We find that in the arm regions, am is flatter than that in the in¬ 
terarm regions. This was already indicated in Basu et al. (2012) and 
Tabatabaei et al. (2007). There appears to be a sharp transition of 
spectral index above and below ~ 50 MqPC^^. We find the spectral 
index to remain constant or get flatter with increasing Sgas above 
~ 50 Mopc^^. There is a sharp steepening of the spectral index be¬ 
low 50 Mopc^^, typically in outer parts and interarm regions. We 
find am close to the injection value of -0.5 for regions having 


Egas in the range 50 - 250 Mopc At even denser regions, we see 
the ttnt to be flatter than -0.5 for the galaxy NGC 5236. 

The error on am depends on the signal-to-noise ratio (S/N) 
of the flux densities within each synthesized beam at each radio 
frequency. The regions of high gas density (Egas > 50 Mopc^^) are 
also typically the regions with strong radio emission with S/N > 5. 
Such regions have < 5 percent error in the a„,. While the low gas 
density regions (Egas < 50 Mopc““) with radio flux density S/N < 5, 
am have up to 20 percent error. Additional error on the estimated 
value of a^t is incurred due to uncertainties in the estimated thermal 
emission. This is primarily due to the unknown value of the electron 
temperate (T,,), assumed to be lO'* K, while estimating the thermal 
emission (see Tabatabaei et al. 2007). As a result, the uncertainty 
in the thermal fraction (fit) is ~ 10 and ~ 15 percent at 0.33 and 
1.4 GHz, respectively (Tabatabaei et al. 2007; Basu et al. 2012). 
This introduces an error of ~ 5 percent on the am. However, note 
that the uncertainty in fn, affects am only towards the inner parts 
of the galaxies where ft is high (Basu et al. 2012). The a„, in the 
outer parts are not affected due to low thermal emission. Overall, 
the estimated am have maximum error of ~ 7 percent towards the 
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inner parts of the galaxies, i.e., regions of higher and up to 
~ 20 percent in the outer parts, i.e., regions of lower Egas. We show 
representative errors on a„i for some of the points covering the span 
of Sgas for each galaxy in Figurej^ 

We checked the possibility of free-free absorption of the syn¬ 
chrotron emission by thermal electrons giving rise to the observed 
flattening at lower frequencies. The observed flux density can be 
expressed as Sy oc v““‘ exp(-Tff). Tg is the free-free optical depth at 
a radio frequency v and is given by. 


Tg = 0.082 



Flere, Tg is the temperature of the thermal electrons assumed to be 
10"* K and EM is the emission measure. The emission measure for 
our sample galaxies are in the range of few 100s (in the disk and 
outer parts of the galaxies) to lO'* cm^^pc (towards the center and 
Fill regions) (Basu et al. 2012). For this range of £Ms, the Tg is ex¬ 
pected to become 1, below 0.1 GFIz. At 0.33 GFlz, Tg for our sam¬ 
ple galaxies lies in the range 10“^ to few times 10“^. Moreover, for 
free-free absorption to affect at 0.33 GFlz, EM x 10® cm“®pc is re¬ 
quired. Regions with such extreme £Afs are rare in galaxies except 
perhaps in the galactic centers (Roy & Rao 2004). The flattening of 
the radio continuum spectra at 0.33 GFlz is therefore unlikely to be 
caused by free-free absorption. 


4 DISCUSSION 


We have studied the variation of the non-thermal spectral index 
with the surface mass density of neutral gas for four spatially re¬ 
solved nearby galaxies. Our results indicate that the local ISM gas 
density plays an important role in shaping up the local radio con¬ 
tinuum spectral behavior. It is therefore difficult to interpret the ob¬ 
served galaxy-integrated power-law nature of the spectrum, as it 
averages the clumpy ISM. As pointed out earlier, after the CREs 
are injected into the ISM, they are subjected to diffusion and en¬ 
ergy loss mechanisms that depend on the in-situ ISM density. De¬ 
pending on the dominant loss process, the injection spectrum of 
the CREs changes differently. Synchrotron and inverse-Compton 
losses have the same effect of steepening the spectrum. The syn¬ 
chrotron timescale (4yn) for CREs emitting at a radio frequency 
Vsyn and having energy E [a; 7.9(vsyn/GHz)'/^(Sj^/yuG)“'^^ GeV] is 
given by. 


= 8.35 X 10® 



yr. 


(5) 


Here, is the magnetic field strength in the plane of the sky. 

The timescale for inverse-Compton losses (?ic) due to scatter¬ 
ing of CREs by radiation field in the galaxy with energy density 
Urad and cosmic microwave background with energy density Ucmb 
is given by. 


he = 4.5 X 10* 


( ^ V 

' / Urad + UcMB \ 

\GeV/ 

\ 10“'2 erg cm“3/ 


yr. 


(6) 


Ionization losses flattens the spectrum towards lower radio fre¬ 
quencies and the timescale (/]„„) is given by (Murphy 2009), 


= 4.1 X 10® 


( 


[3 lnf-^)-t42.5l 

\ cm“3 j 

\GeV/ 

[ \GeV/ J 


yr. (7) 


Here, (n) is the average number density of neutral gas (atomic -l- 
molecular) as seen by the CREs. 


Bremsstrahlung losses does not significantly change the spec¬ 
trum and its timescale (tbrem) is given by (Murphy 2009), 

rtam = 3.96xl0^f-^) yr. (8) 

\cm 

Since the ISM of a galaxy is not smooth in its physical prop¬ 
erties, i.e., magnetic field strength and gas density, different losses 
will dominate in different parts of the galaxy. To compare the var¬ 
ious energy loss scenarios we plot these energy loss timescales as 
a function of Egas in Figure [3 As these loss timescales depend on 
the energy of the CRE and thereby the radio frequency, we show 
these timescales at 0.1, 0.3, 1.4, 5 and 10 GHz. Further, at a fixed 
frequency, the effective CR energy depends on the magnetic field. 
The top axis in Figure [3] shows the energy of the CREs. The mag¬ 
netic field strength (S) can depend on Egas as B oc E’aj (Schle¬ 
icher & Beck 2013). Using the normalizatior[^ given by Schleicher 
& Beck (2013), B x 2 x (Egas/Mgpe^^)'/^ fiG. Further, for the 
inverse-Compton losses, we considered the t/„d to depend on Egas 
as, f7rad = 8.33x 10“'®(Egas/Mopc“^) erg cm“*, normalized in terms 
of the Milky Way (Schleicher & Beck 2013). 

The average number density of total neutral 
(atomic-l-molecular) gas ((«)) was estimated from surface 
mass density (EgaJ. Typically, the surface mass density of atomic 
gas as traced by Hi emission does not exceed ~ 10 Mopc^^ 
for normal star-forming galaxies (see e.g., Leroy et al. 2008; 
Kalberla & Dedes 2008). We therefore assumed the regions with 
2igas <: 25 Mopc^^ to be originating entirely from molecular gas 
as traced by CO. Such regions are typically found in the inner 
regions of the galaxies and in spiral arms. The regions with 
2ga.s £ 10 Mopc^^, were assumed to be originating from atomic 
Hi gas, typically found in interarm regions and outer parts of the 
galaxies. For regions having Egas ^ 10-25 Mopc^^, we assumed a 
mixture of atomic and molecular gas such that Hi/(Hn-H 2 ) varying 
linearly between 1 and 0. 

The surface gas density is converted to mid-plane number den¬ 
sity assuming a Milky Way type model for the vertical distribution 
of the atomic and molecular gas. We used a 3-component vertical 
distribution for the Hi, following Dickey & Lockman (1990) and 
a single component for H 2 , following Sanders, Solomon & Scov- 
ille (1984). The three Hi components are as follows: 1) a Gaussian 
component with FWHM of 212 pc (cri = 90.2 pc) having mid¬ 
plane density of ni, 2) a second Gaussian component with FWHM 
of 530 pc ((T 2 = 225.5 pc) having mid-plane density of 112 and 3) 
an exponential component with scale height (1) of 403 pc with mid¬ 
plane density (Dickey & Lockman 1990). The components have 
relative ratios, ni : 112 '■ 112 = 0.395 : 0.107 : 0.064 (Dickey & 
Lockman 1990). Thus, Ehi is modeled as. 


+00 

Jim = J' [«i exp(-z^/2a-j) - 1-112 exp(-z^/ 2 a- 2 ) 

+ n3exp(-z//)]dz (9) 

Similarly, the mid-plane number density of H 2 («h 2 ) is com¬ 
puted assuming a FWHM of 110 pc (Sanders et al. 1984) corre- 


^ Note that the normalizations in Schleicher & Beck (2013) are given in 
terms of the surface star-formation rate (Ssfr)- We converted SsFR to Sgas 
following Kennicutt (1998). 
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Figure 3. The various CRE energy loss timescales as a function of the neutral gas surface density (Xgas). The different panels are for CREs emitting at 0.1, 0.3, 
1.4, 5 and 10 GHz. The shaded region shows the range in Sgas where the spectral index is observed to change (see Figure[2j- For comparing the CRE injection 
timescales, we plot the supernova surface rate density (number of supernova per pc“) as dotted lines. 
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spending to cthj = 46.8 pc. Thus, Eh, is given by, 

+00 

Sh, = exp(-z^/2crH,) rfz (10) 

— CO 

The total mid-plane number density («tot) is given by, ntot = 
ni + ^2 + ”3 + uh,- We assumed the average number density of 
neutral gas, (n), seen by the CREs to be half of «tot- 

4.1 CR injection 

All the CR energy loss mechanisms would be relevant if the CR 
injection timescales were significantly larger. Otherwise, the loss 
mechanisms may not affect the radio continuum spectrum and one 
would observe the spectral index close to the injection value, ai„j. 
Under DSA, depends on the Mach number (M) as, Ofjnj = 
-(M^+3)/(2M^-2) (Blandford & Eichler 1987). For strong shocks, 
ffinj ~ -0.5. We note that the Mach number depends on the sound 
speed and thereby the gas density (p) as M oc yfp. Thus, depending 
on the gas density, ttinj could be in the range-0.5 to -0.8. 

The black dotted lines in Figure [3 show the injection 
timescales of freshly generated CR particles estimated from the 
supernova rates. The stars that are > 8 Mq end up in Type II and 
Type Ib supernova, and are responsible for acceleration of the CREs 
in normal galaxies (Condon 1992). Using the initial mass func¬ 
tion (IMF) given by Kroupa (2001), we estimated the mass frac¬ 
tion (/m) of stars that could form supernovae by integrating the 
IMF between 8 and 150 Mo. Thus, for a surface mass star forma¬ 
tion rate of Esfr, the mass rate of forming supernova stars is given 
/m X 2sfr. The mean surface supernova rate (vsn) is given by, 
''SN = /m X Esfr/(Msn). Here, (Msn) is the mean supernova mass. 
In our case, /m/(Msn) is found to be 0.01. Esfr was estimated from 
Egas using the Kennicutt-Schmidt law (Kennicutt 1998). From Fig- 
ure[3we find the local CR injection timescales to be larger than the 
loss timescales and hence, the freshly generated CREs do not con¬ 
tribute significantly to the observed synchrotron spectrum. Except 
for regions of high star-formation or recent starburst activity, an ob¬ 
served flatter spectral index of -0.5, close to the injection value, is 
unlikely to be caused by the fresh CREs. 

4.2 Gas density and energy losses 

From Figure it is clear that ionization and bremsstrahlung losses 
dominate in regions of high gas density at lower radio frequencies 
(< 0.3 GHz) where the synchrotron emitting CREs have less en¬ 
ergy (~ 1 GeV). Inverse-Compton losses are mostly un-important 
within galaxies, except for very low density regions with magnetic 
field strength < 3.3 pG. For regions with Egas > 50 Mgpe^^, hon 
and tbiem are lower than that of below ~ 0.3 GHz. Hence, such 
regions are expected to have flatter spectral index as is observed 
in Figure |2] (shown as shaded area in Figure [Sj. The regions with 
Sgas ^ 200 MqPC^^ are perhaps contributed by dense molecular 
clouds having high densities and magnetic field strengths. In such 
regions, ionization losses can become significant at 0.3 GHz, if 
the magnetic field strengths are high, giving rise to spectral index 
flatter than their injection value of ~ -0.5. This is seen for the 
galaxy NGC 5236, where regions with Egas <: 200 Mgpe^^ are ob¬ 
served to have in the range -0.3 to -0.5. Such regions have been 
observed to have relatively high magnetic field strengths (> 25 pG; 
Basu & Roy 2013). Bremsstrahlung losses may, however, dominate 
significantly in moderately dense regions with >100 Mgpe^^ at 1.4 
GHz. Such regions do not affect and hence, we do not expect 


to observe flatter than -0.5 when measured between radio fre¬ 
quencies above 1.4 GHz. For ionization losses to dominate over 
synchrotron losses at higher radio frequencies (between ~ 1.4-5 
GHz), an extremely high surface gas density of > 1000 Mopc^^ is 
required. However, in such regions, the injection timescale of CREs 
is comparable to the loss timescales within a factor of 2. 

Regions with Eg^s 5 50 Mgpe^^, are observed to have steep a„, 
(see Figurel^. From Figure^ we find that such regions are dom¬ 
inated by synchrotron and/or inverse-Compton losses at 1.4 GHz 
and hence perhaps causes steepening of «„!• However, note that the 
CRE injection timescale approaches few times 10® yr in regions of 
even lower Eg^^ (< 10 Mepc^^), typically found in the interarm re¬ 
gions and outer parts of the galaxies. The CRE population in such 
regions are dominated by propagation from star-forming regions in 
adjacent arms. Depending on the propagation mechanism, i.e., sim¬ 
ple diffusion or streaming instability at Alfven velocity, the CREs 
at 0.33 GHz propagate 1.4 to 2 times longer distances than those 
at 1.4 GHz (see Basu & Roy 2013). This could also give rise to 
the observed steepening of a„i between the two frequencies in such 
regions. 

To test the scenario of the expected differences in spatially re¬ 
solved ttnii we studied between a low frequency pair (0.3 and 
1.4 GHz) and between a high frequency pair (1.4 and 4.85 GHz). 
The galaxy NGC 6946 has the best resolution spectral index map in 
our sample. NGC 6946 was observed at 4.85 GHz using the VLA in 
D-array configuration and a combination with a single dish Effels- 
berg observation was done to ensure no missing flux density (Beck 
2007). In Figure|4](left-hand panel) we plot the total radio emission 
spectral index between 1.4 and 4.85 GHz [a(1.4-4.85 GHz)] versus 
the spectral index between 0.33 and 1.4 GHz [Qr(0.33-1.4 GHz)]. 
The total emission spectral index estimated between the two sets 
of well separated radio frequencies follow the 1:1 line closely. The 
symbols are colour coded based on logj,, Eg^j. In the high gas den¬ 
sity regions, typically in the arms and inner regions of the galaxy, 
the spectral indices are flatter. However, these values could be af¬ 
fected due to the presence of high thermal emission. 

In Figure |4] (right-hand panel), we plot the non-thermal spec¬ 
tral index estimated between 1.4and4.85 GHz [ant(1.4-4.85 GHz)] 
and between 0.33 and 1.4 GHz [a„,(0.33-1.4 GHz)] (see Section 2). 
In general, Q'nt(l.4-4.85 GHz) is found to be steeper than Q'nt(0.33- 
1.4 GHz). It is steeper by up to -0.5 in regions of lower gas surface 
density (bluer points), i.e., the interarm regions and outer parts of 
the galaxy. This is consistent with the scenario that such regions are 
dominated by synchrotron and/or inverse-Compton losses resulting 
in steeper spectral indices at higher frequencies. 

We compare the expected variation of a„i measured between 
the two pairs of radio frequencies using the different models of en¬ 
ergy loss, i.e., the continuous injection (Cl), Jaffe-Perola (JP) and 
Kardashev-Pacholczyk (KP) models. We generated synthetic radio 
spectra for each of the models with varying break frequencies. In 
Figure |4] (right-hand panel) we plot the trajectory of for a con¬ 
stant injection spectral index of aj„j. The Cl model is indicated by 
the dashed-dot lines, the JP model is indicated by dashed lines and 
the KP model is shown as the dotted lines. It is clear that the Cl 
model cannot give rise to the observed steepening of the spectrum 
between the two pairs of radio frequencies of our study. This is also 
indicated in Figure[3] which shows that the injection timescales are 
larger than the loss timescales. The JP and KP models reproduce 
the observed steepening at higher frequency better. The maximum 
steeping observed in our case corresponds to a break frequency of 
~ 4 GHz for the JP model and ~ 2.5 GHz for the KP model. De¬ 
tailed modelling of the spatially resolved radio continuum spectra 
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Figure 4. Left-hand panel: Total spectral index between 1.4 and 4.85 GHz 4.85 GHz)] versus spectral index between 0.33 and 1.4 GHz [Q'total(0-33- 

1.4 GHz)] for the galaxy NGC 6946. Right-hand panel: Non-thermal spectral index between 1.4 and 4.85 GHz [Q;nt(l-4^.85 GHz)] versus between 0.33 
and 1.4 GHz [Q'nt(0.33-1.4 GHz)]. The spectral indices are measured within one synthesized beam. The star shows the galaxy integrated spectral index. The 
symbols are colour coded based on the logjg Sg^s. The dashed, dashed-dot and the dotted lines shows the Cl, JP and KP models, respectively, for different 
of -0.5, -0.7 and -0.8 (starting from top right). The solid lines shows the 1:1 line. The symbols with errors represents the typical error on the spectral index 
for different values of log|Q Sgas. 


would require high quality maps at several frequencies covering a 
wide range. 

Moreover, we note that the regions of high gas density are also 
the regions with dominant ionization and/or bremsstrahlung losses. 
Such losses have the effect of flattening ctni at lower frequencies as 
seen in Figure [2] As is also evident from Figure [3] at regions of 
high surface gas density (> 200 Mopc^^), ionization losses should 
start dominating at ~ 0.3 GHz. Such regions are observed in the 
galaxy NGC 5236 and are likely to be giant molecular clouds. We 
observe ctni to be flatter than the typical injection values of -0.5. 
However, these losses will not affect at higher radio frequencies. 
We therefore observe the values of Qr„t(l.4^.85 GHz) to be well 
mixed in all regions with no clear trenqj with gas density. 

Due to the smooth transition of the values of ttni at lower and 
higher frequencies in the spatially resolved case, the galaxy inte¬ 
grated case does not capture this scenario well, i.e., no significant 
change in spectral index should be expected. The stars in both the 
panels in Figure |4] shows the galaxy integrated spectral index. The 
spectral index for total radio emission is observed to be close to 
the 1:1 line. For the non-thermal emission (right-hand panel), the 
integrated spectrum is still observed to remain close to the 1:1 line. 
Within errors, we observe both total and non-thermal spectral index 

to remain close to-0.8 for both pairs of frequencies. Hence, we 

do not see any significant change of the slope for the galaxy inte¬ 
grated spectrum. Moreover, the dense regions that shows evidence 


'* Note the systematic trend of density vaiiation (red to orange to yellow to 
blue points) as one moves from right to left along x-axis. No such trend is 
observed while moving from top to bottom along y-axis. 


of ionization losses, cover only a small fraction (< 5 percent) of the 
galaxy. Thus, the galaxy integrated spectrum may not show such 
a behavior. Our study comprehensively shows that different loss 
mechanisms are dominant in different parts of the galaxies. It is 
difficult to single out any one of the energy loss mechanisms while 
interpreting the galaxy integrated spectrum. Such studies, hence do 
not often show significant deviation from a power-law behavior for 
a wide range of radio frequencies. 

Data on more galaxies are required to see such clear trends. 
LOFAR and EVLA are ideal instruments to do such studies ro¬ 
bustly at lower (< 150 MHz) and higher (> 1.4 GHz) radio frequen¬ 
cies, respectively. Complementary data from the GMRT covering 
the intermediate frequency range (150 - 610 MHz) could help us 
model the spatially resolved spectral behavior of continuum emis¬ 
sion in detail. High resolution CO observations from the ALMA 
could facilitate measuring the distribution of the molecular gas at 
unprecedented detail. 


5 CONCLUSIONS 

We have studied the variation of spatially resolved non-thermal 
spectral index (arm) as a function of the total neutral gas surface 
mass density (Sg^O for four nearby galaxies. The galaxies in our 
sample are, NGC 4736, NGC 5055, NGC 5236 and NGC 6946. 
In this study, we have measured the and Zg^s at sub-kpc linear 
scales except for the galaxy NGC 5236. The spectral index were 
measured between 0.33 and 1.4 GHz. We present our conclusions 
here. 

(i) The local ISM gas density plays an important role in shaping 
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up the radio continuum spectrum. We find the to remain 
comparatively flat (> -0.7) in the regions of high gas density 
(Sgas <: 50 Mopc^^) and steepens sharply (< -1) in regions of 
low gas density (Sg^s < 50 MqPC^^). 

(ii) The observed flattening is unlikely to be caused by free- 
free absorption at low frequency such as 0.33 GHz. We 
find that in the regions of high gas density (> 50 Mepc^^), 
bremsstrahlung and ionization losses could flatten the a^i at 
radio frequencies below ~ 0.33 GHz. 

(hi) Synchrotron losses dominates at higher radio frequencies 
(> 5 GHz) in regions having moderate densities (Sgas < 
100 MqPC^^). In denser regions, bremsstrahlung losses are ex¬ 
pected to dominate. Inverse-Compton losses are mostly unim¬ 
portant in galaxies except for regions of very low gas densities 
having magnetic field strengths < 3.3 /rG. 

(iv) For the galaxy NGC 6946, we study the variation of q:„, mea¬ 
sured between 0.33 and 1.4 GHz and between 1.4 and 4.85 
GHz. We find that between the higher frequency pair is 
steeper than that between the lower frequency pair indicating 
synchrotron or inverse-Compton losses at 4.85 GHz. The ob¬ 
served steepening cannot be caused by continuous injection 
of CR particles and requires single shot particle injection. 

(v) Due to the clumpy nature of the ISM, these local CRE en¬ 
ergy losses may not show indications in the overall galaxy- 
integrated spectrum. Thus, the galaxy-integrated spectrum 
may remain a power-law for a wide range of radio frequen¬ 
cies. 
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